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Abstract. The canonical quantization of macroscopic electromagnetism was recently 
presented in New J. Phys. 12 (2010) 123008. This theory is here used to derive 
the Casimir effect, by considering the special case of thermal and zero-point fields. 
The stress-energy- momentum tensor of the canonical theory follows from Noether's 
P5 , theorem, and its electromagnetic part in thermal equilibrium gives the Casimir 

^ ' energy density and stress tensor. The results hold for arbitrary inhomogeneous 

^H. magnetodielectrics and are obtained from a rigorous quantization of electromagnetism 

in dispersive, dissipative media. Continuing doubts about the status of the standard 
Lifshitz theory as a proper quantum treatment of Casimir forces do not apply to the 
^ , derivation given here. Moreover, the correct expressions for the Casimir energy density 

' and stress tensor inside media follow automatically from the simple restriction to 

thermal equilibrium, without the need for complicated thermodynamical or mechanical 
arguments. 



PACS numbers: 42.50.Lc, 42.50.Nn, 12.20.-m 
1. Introduction 

The phenomenon of forces on macroscopic objects due to electromagnetic zero-point 
and thermal fields is usually called the Casimir effect, after the first, highly idealized, 
theoretical result in the subject A theory of the effect for realistic materials was 
given by Lifshitz and co-workers [21 [3l H] and this remains the most general formalism for 
describing the local electromagnetic quantities that produce the forces. Some authors 
prefer other terminology, such as Casimir-Lifshitz effect, or van der Waals forces, to 
designate this same phenomenon; in this paper the term Casimir effect is used, both for 
the zero-point and thermal contributions. 

Despite its good agreement with experiments that measure the zero-point Casimir 
force, there have been been persistent doubts about the status of Lifshitz theory as a 
quantum theory. These doubts are voiced, for example, in two recent publications [5l|6], 
and a detailed analysis in [6] concludes that "Lifshitz theory is actually a classical 
stochastic electrodynamical theory" . The source of these doubts is not difficult to find 
if one examines the details of Lifshitz theory. As the Casimir effect is a phenomenon of 
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quantum electromagnetism in the presence of macroscopic media, one would expect the 
general theory of the effect to be based on the principles of quantum electrodynamics. 
Lifshitz theory, however, is based rather on the the principles of thermodynamics. In 
fact there is no Hamiltonian in Lifshitz theory and there are no quantized fields, yet the 
formalism seeks to calculate the forces caused by quantum zero-point fields (as well as 
by thermal fields). 

The requirements of a quantum theory of light in media that can serve as a basis for 
the Casimir effect (among many other phenomena) are easy to identify. The theory must 
hold for arbitrary magnetodielectrics in order to be applicable to the realistic materials 
used in experiments, materials whose optical properties (dielectric functions) are known 
only through measurement. As the Casimir effect is a broadband phenomenon in which 
all frequencies must be included, the theory must also take full account of material 
dispersion and absorption. The classical theory of light in these circumstances is of 
course the macroscopic Maxwell equations, where the electromagnetic properties of the 
media are encompassed in arbitrary electric permittivities and magnetic permeabilities 
obeying the Kramers-Kronig relations. The theoretical basis of the Casimir effect should 
therefore be a quantum theory of Maxwell's macroscopic electromagnetism. But whereas 
the quantum theory of the free-space Maxwell equations, quantum electrodynamics 
(QED), was the foundation of quantum field theory, and led to a general formalism 
for quantizing classical field theories, this quantization procedure was not applied 
to macroscopic electromagnetism. What has become known in the literature as 
macroscopic QED is not based on the rules for quantizing field theories, but is instead 
a phenomenological theory wherein no rigorous quantization is attempted (see [71 [8] 
for detailed presentations). This phenomenological procedure is subject to much of 
the criticism directed at Lifshitz theory, and it will be seen that the results derived 
in this paper require, among other things, an action principle, something that is 
lacking in the phenomenological approach. The quantization rules of quantum field 
theory were not obviously applicable in the case of macroscopic electromagnetism 
due to the complications of dispersion and dissipation, and this is what led to the 
phenomenological approach. In fact it had been "widely agreed" [9] that a proper 
quantization of macroscopic electromagnetism could not be performed, and that only 
in cases where a simple microscopic model of the dielectric functions of the medium is 
explicitly introduced could the standard quantization rules be applied. In [10], however, 
the canonical quantization of macroscopic electromagnetism was achieved, providing a 
rigorous macroscopic QED and removing the need for a phenomenological approach. 
Since the macroscopic QED derived in [10] applies to arbitrary magnetodielectrics and 
takes full account of dispersion and absorption, it meets the criteria outlined above for 
a rigorous quantum foundation for the Casimir effect. This paper derives the Casimir 
effect from macroscopic QED by considering the special case of thermal equilibriumj^l 

A brief summary of macroscopic QED is given in section [2l The stress-energy- 

:j: Hereafter macroscopic QED refers to the canonically quantized theory presented in [10 , not the 
previous phenomenological formalism (Tj |8] . 
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momentum tensor of macroscopic electromagnetism is derived in full generality in 
section [3|, by application of Noether's theorem to the action principle given in [TO] . 
Specialization to thermal equilibrium (including the zero-point fields) is made in 
section HI The expectation value of the electromagnetic part of the stress-energy- 
momentum tensor in thermal equilibrium gives the Casimir effect. Correlation functions 
of the quantum field operators in thermal equilibrium are calculated in [5] and used in 
sections [6] and [7] to obtain the Casimir energy density and stress tensor. 



2. Macroscopic QED 



This section summarizes the results of macroscopic QED that we require to derive 
the Casimir effect. The action of macroscopic electromagnetism is [10] 

A, X^, YJ = S,r^[cP, A] + 5x[X^] + 5y[Y^] + 5int[0, A, X^, Y^], 

where S'em is the free electromagnetic action 



5'em[0, A] 



^ / d^x -E ■ E - B ■ B , Ko = 



(1) 

(2) 



and Sy are the actions for free reservoir oscillators: 



^y[Y^] = lj^^^l da; {dtY^ ■ dtY^ - co^Y^ ■ Y^) , 



(3) 



(4) 



and ^int is the interaction part of the action, coupling the electromagnetic fields to the 
reservoir: 

5i„t[0, A, X^, Y^] = y d's; ^ da; [a{r, ■ E + /3(r, a;)Y^ ■ B] , (5) 

2ko 



a r, a; 



2e 







.1/2 



vr 



-a;£:il^r,a;j 



/3(r,a;) 



vr 



1/2 



(6) 



The imaginary parts of the dielectric functions of the medium appear in the coupling 
functions (Q; their real parts are given by the Kramers-Kronig relation [TTl [12] : 



£R(r,a;')-l = -P / 

71" Jo 



^2 _ 



and similarly for K(r,u). (7) 



As in [TU], we assume the medium is isotropic, with scalar dielectric functions e{r,u) 
and K,(r,u) = l//i(r,a;); anisotropy can be included by obvious modifications. The field 
equations of the action ([I])-(l6]) are 

POO 

£oV-E+ / da;V- [a(r,a;)Xj = 0, (8) 
Jo 

POO 

-KoV xB + eodt'E+ da; {a(r, a;)9tX^ + V x [/3(r, a;)Yj} = 0, (9) 

Jo 



'0 

a(r,u)E = 0, 



-d^Y^ - u^Y^ + (3{r,u)B = 0. 



(10) 
(11) 
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By solving the equations for the reservoir fields and Y^^, the independent equations 

of the electromagnetic fields are found and are precisely the macroscopic Maxwell 
equations [10] 

V ■ D = a, (12) 

VxH-9iD=j, (13) 

where the charge density a and current density j are given in terms of arbitrary free- 
field solutions for the reservoir fields X^^ and Y^^. The other two Maxwell equations are 
identities in terms of the potentials A and 0. 

As the action ([I])-(l6]) features the dynamical fields and their first derivatives, 
canonical quantization can proceed without difficulty [in]. The resulting Hamiltonian 
can be diagonalized to the form 



H 



A=e,m 



d'r 



dw f\wC\{Y, u) ■ C'x{r, u) 



(14) 



0.(15) 



where the diagonalizing eigenmode creation and annihilation operators obey 

Cxi{r,u),Cl,-{Y\uj') = Sij6x\'S{u - u')S{r - r'), Cxi{r,uj),Cx'j{r' 

The infinite zero-point energy has been omitted from the Hamiltonian ( fT4l) . Zero-point 
energy is of course crucial for the Casimir effect, but the electromagnetic energy density 
and stress tensor will be calculated below taking full account of the zero-point fields. 
The quantum macroscopic Maxwell equations 

V-D = a, (16) 

VxH-9iD=j, (17) 

hold with the charge and current density operators given, in the frequency domain, by 



(3"(r, oj) = —V ■ j(r, u) = — 27rV 



Jir,a;j 



-27ria; 



he, 







nl/2 



TT 



-ei{r,u 



CJr,uj) + 27rV x 




nl/2 



Cm(r,w) 



.(19) 



The relationship between fields in the time and frequency domains is, for the example 
of the electric field, 

r"00 



E(r,t) 



1 
2^ 



du 







E(r, u) exp(— iwt) -|- c.c. 



(20) 



The electromagnetic field operators and the reservoir field operators are also expressed 
in terms of the diagonalizing operators, as follows. For the electric field operator the 
relation follows from ( !T9|) and 



E(r,t) 



2tt 



du I dh' 



ici;G(r, r', u) ■ j(r', u) exp(— iwt) -|- h.c. 



where the Green bi-tensor G(r, r', w) is the solution of 



OJ 



V X [K{r, u)V X G(r, r , u)] - —e{r, w)G(r, r', u) = 16{r - r'). 



(21) 



(22) 
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The magnetic field operator is 

B(r, u) = -iV X E(r, u) /u. (23) 

Finally, the reservoir field operators and can be written in the frequency domain 

as 



X^(r, cj') = 2tt 



h 




—6{io-uj')C,{r,uj') 



Y^(r,a;' 



+ 



2n 



a r, u 



2uj 

h 



+ 



1 




2u 
i/3(r,u; 
2uj'uj 



iO+ u + oj' 
5{uj-uj')C,,,{r,uj') 



E(r,u;'), 



(24) 



1 



1 



V X E(r,a;') 



(25) 



w — uj' — iO+ u + oj' 
The last two equations were not explicitly written in [TOj; they are obtained from the 
expansion (58) in [TO] of the operator X^ in terms of the diagonalizing operators, and 
the analogous expansion of Y^, by inserting the results (70), (71), (76) and (88) of p^ . 
Use is also made of the fact that all frequency arguments are taken only at positive 
values (see (l20l) ). so that the sums of frequencies in denominators do not give rise to 
poles. 



3. Stress-energy-momentum tensor of macroscopic electromagnetism 

The canonical formulation of macroscopic electromagnetism given in [10i| leads directly 
to the stress-energy-momentum tensor of the system, in both the classical and quantum 
cases. As dissipation of electromagnetic energy by the medium is fully taken into account 
through the presence of the reservoir fields, a conserved total energy of the system 
exists (if the dielectric functions of the medium are time independent). Similarly, total 
momentum is conserved if the medium is homogeneous. 

The complete stress-energy-momentum tensor of the system follows from 
application of Noether's theorem to the action ([I])-(l6l). We state the results for classical 
fields, but the quantum stress-energy-momentum tensor is of the same form because the 
quantum field operators obey the same dynamical equations as the classical fields, as 
shown in [TO] . 

3.1. Energy density and energy flux 

The invariance of the action ([I])-(l6]) under active time translations of the dynamical 
fields implies a conservation law that can be extracted as follows [13J. We make an 
active infinitesimal time translation of all the dynamical fields — in the case of the vector 
potential, A(r, t) — )■ A(r,t + (^(r, t)) — but take the translation C(r, i) to vary in space 
and time. The resulting change in the action can be reduced to the form 

5S= f d^x{pdtC + s-WC), (26) 
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where p is the energy density and s is the energy flux, obeying the conservation law 

+ V ■ s = 0. (27) 
The calculation is straightforward and yields 



p= y 



+ 



3E ■ {-dtA + V0) + B 
dcj 



c 
00 







aX^ ■ V0 - /3Y^ ■ B 



1 1 /"°° 

-E'dt<j) + 2V^'A^^dtA, - / da; [aXldt(j) - /SY^e/'dtAk] 
c _ Jo 



f28) 



(29) 



where anti-symmetrization of tensor indices is denoted by square brackets and e^^^ is 
the Levi-Civita tensor [14j. As in the case of free-space electromagnetism (121 |T3], the 
energy density (!28l) and flux (129|) that directly emerge from Noether's theorem are not 
gauge invariant. They are however equivalent to gauge-invariant quantities because they 
fail to be gauge invariant up to terms that identically satisfy the conservation law ([2] 
Specifically, the quantities 



ri t 
J t 



t i 



(30) 
(31) 



-2ko0VM^i + 
identically satisfy 

^N^f\' + V.{dtf: + V,f7) = 0. (32) 

Comparing fl32|) with fl27|) . we see that if Vif\^ is added to p, and (9t/*/ + V jf\^ is 
added to s*, then the conservation law ( 1271) will still hold. With use of the field equations 
(!8|)- (fTT]) . the energy density and flux that result from these additions are gauge invariant 
and are given by 

1 



P 



2 

+ 



-E^ + B^ 



00 







s = kqE X B 



dul3E X Y,^. 



Yl)-/3Y^-B 



(33) 
(34) 



It is straightforward to verify that the conservation law f l27l) holds for (1331) and 
when the fields obey the dynamical equations fl8l)- f|TTl) . 



3.2. Momentum density and stress tensor 

The total momentum of the electromagnetic field plus matter is conserved in flat space- 
time. In the description of macroscopic electromagnetism, however, the microscopic 
degrees of freedom of the magnetodielectric medium are not included in the action; as 
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a result, the translation symmetry that gives rise to momentum conservation is not in 
general present in the dynamical system considered. A conservation law for momentum 
will exist within the macroscopic framework only if the coupling functions ([6]) are 
independent of position (a homogeneous medium) so that the action is invariant under 
active spatial translations of the dynamical fields. It is nevertheless instructive to retain 
the general case of an inhomogeneous medium; the conservation law for momentum 
will then be seen to fail in the inhomogeneous case due to the appearance of spatial 
derivatives of the coupling functions (E]). 

We make an active infinitesimal spatial translation of all the dynamical fields; in 
the case of the vector potential this is A(r, t) — t- A(r + w(r, t), t). The resulting change 
in the action can be written 

(55 = - y d^x {pi dtw' + a/ Vjw') , (35) 

where the momentum density p and stress tensor a^-' obey, in homogeneous media, the 
conservation law 

dtPi + ^jcr^' = 0. (36) 
The form (l35l) is achieved with 

POO 

SoE^W.A, - / du {dtXiW.X^j + dtYiViY^, - aXiV^A,) , (37) 

C5/ + /to [^E^^^'P + 2V[M'=lV.Afc^ ^ / "^"^ («^iV.0 - PY^eJiV^A') , (38) 

where C is the Lagrangian density, i.e. the integrand in the action ([1]). Again, the initial 
results ( 1371) and ( l38l) are not gauge invariant. The quantities 



p. 



/V := -e,A,E^ - / da; aA.Xl =: -/*/, (39) 

POD 

/V := -2ft:oAiV[M'=l + / duPe'^'A^Y^i, (40) 
identically satisfy 

WV + V,WV + V./V) = 0. (41) 

Thus, addition of Vjf''/ to pi, and of i9t/*/ + Vkf^/ to a/, does not affect the 
momentum conservation law fl36|) . After these additions, and use of the field equations 
(l8|)-(fTT1). we obtain a gauge-invariant momentum density and stress tensor: 

POO 

Pi = £o(E X B), - / da; [dtXlV^X^, + 9tFj + a(B x X^),] , (42) 
^^ = ^-^/'(^oE^ + «:oB2) - e,E,E^ - t^oB.B^ 

+ / do; <! 6/ -{d,X^f + -{dtY^r - -uj'i^l + Yl) + aX^ ■ E 

-aE,Xi + pY^,B^]. (43) 
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When the field equations (15])- (ITT]) hold, the momentum density (1421) and stress tensor 
fH3|l satisfy 



so that the conservation law (!36l) indeed holds for homogeneous media. 
4. Thermal equilibrium 

Casimir forces can be calculated from either the electromagnetic energy density or stress 
tensor in the presence of macroscopic media. General electromagnetic fields will exert 
forces on the media, but the Casimir effect is the special case when the fields are in 
their ground state (zero-point fields) or in thermal equilibrium with the media (the 
latter case of course includes the contribution of the former). To derive the Casimir 
effect we therefore need only assume that the eigenmodes of macroscopic QED are in a 
thermal mixed quantum state. The expressions for the electromagnetic energy density 
and stress tensor that determine the forces then follow from the general results fl33|) and 
( 143|) (which, as noted above, also hold for quantum field operators). This is in sharp 
contrast to Lifshitz theory [21 [3l H] , where extraordinarily complicated thermodynamical 
and mechanical arguments are required to obtain the electromagnetic stress tensor in 
media and vacuum, in a manner that has nothing obvious to do with the principles 
of quantum mechanics [6]. Although the calculations below are certainly tedious, the 
only ingredients are the quantum theory of macroscopic electromagnetism [lOj and a 
restriction to the case of thermal equilibrium. 

To impose thermal equilibrium on the bosonic eigenmodes of macroscopic QED, we 
assume that the expectation value of the number operator of the eigenmodes is given 



Each eigenmode at each frequency oo and each position r is a quantum harmonic 
oscillator and so the expectation value for the number of quanta (excitation level) 
in each of these oscillators should be given by the Planck distribution (1471) . The 
complication in implementing this simple prescription is that the eigenmodes are a 
continuum in frequency and position. The obvious way of handling this last fact is to 
use delta functions as in f H5|) . a procedure also followed in where a simple quantized 
microscopic model of a medium was analysed. But it must be admitted that there is 
no clear mathematical basis for the density operator which is supposed to underlie the 
expectation value (H5i) . The problem is that the density operator should be defined in 




(44) 



by 
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terms of number states of the quanta, the excitation levels of the oscillators, but because 
a continuum of oscillators is excited, there is no clear way of writing and normalizing 
these number states from which one would then construct the density operator of a 
thermal mixed state. This may well be an inessential technical issue, but one should 
bear in mind that the apparently clear intuition employed in writing (H5ll -( H8ll conceals 
formidable mathematical difficulties. Moreover, careful study of the derivations that 
follow will show that the delta function in frequency in the correlation functions fH5l) 
and ( H6|) is to be regarded as a limit to be taken only in the final stages of the calculations; 
this allows any ambiguities arising from products of delta functions to be negotiated. 



5. Thermal field correlation functions 



As all the dynamical field operators of macroscopic QED are expressible in terms of 
the eigenmode creation and annihilation operators, (I45l)-(l48l) immediately give the 
expectation values of products of field operators (correlation functions) in thermal 
equilibrium. For the current density operator ( IT9l) in the frequency domain we obtain 



jT(r, u) ® j(r', u') ) = AnhAf{u)6{u - u') < u^eoei{r, uj)16{r 



+ KqV X 



— Ki(r, uj) 16{r - r')\/-Ki(r',w')] x V'j (49) 

Ar(.) + l\J>'^'^^J^'^^'^^) ^''^ 
j(r,a;)®j(r',a;')> =0, (51) 



Af{ 



where the notation x V' denotes a curl with respect to the right-hand index, so that 
V(r) X V= V X V(r) for a vector V(r) (note that there is no minus sign included in this 
definition of V acting on the right-hand side).. Equations fH9l) - fl5Tl) can be viewed as an 
example of the fiuctuation-dissipation theorem, but here they are a simple consequence 
of macroscopic QED in thermal equilibrium. 

From (j20l) . the equal-time correlation function for the electric field operator is 
expressed in terms of the frequency- domain correlation function by 

^E(r, t) ® E(r', t)^ = — du du' |exp[-i(a; - u')t] (^E(r, u) ® Et(r', u;')) 

+ exp[i(u; - uj')t] (t^r, to) ® E(r', w')) } • (52) 

(Terms of the form {E{uj)'E{u')) and {W{u)W{u')) vanish.) The frequency- domain 
correlation functions are written in terms of the Green bi-tensor using (120|) . (12T|) and 
(149|) - (!5T1) : after spatial integrations by parts we obtain 

= Anhf^o J dh" ^^icu)u;^5{co - to') |^^i(r", co)GUr, r", u;)G*/(r', r", to) 
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~K,{v\uj)[G{vy,uj)x V'U [G*(r',r",u;)x V"]/| . (53) 

The correlation function {Ei{r,u)Ej{r' ,u')) is given by the complex conjugate of (l53il 
with M{ijj) replaced by Af{uj) + 1. The right-hand side of f l53|) can be simplified as 
follows. We first note the symmetry property of the Green bi-tensor that holds for 
media invariant under time reversal (non-magnetic j§| non-moving media) [1]: 

G,,(r,r',u;) = G',,(r',r,u;). (54) 

Take the matrix product of G*(r", r, u) with ( l22l) and integrate over r; after integration 
by parts and use of fl5^ . the imaginary part of the resulting relation simplifies to 



r'\r',u). (55) 



j dh |-Ki(r, CO) [G*(r", r, u) x VU [V x G(r, r', a;)]^ 

+^£i(r,a;)G*,,(r",r,a;)G^(r,r',a;)| = G 

Use of fl55l) and its complex conjugate, together with fl5^ . allows fl53|) to be simplified 
to 



E^(r, u) E(r', u') ) = Anhfioditu - u')Af{u)u^Gi{r, r', u) (56) 

■^^"^^ ^E{r,u)^E\r',u')). (57) 



A/'(a;) + 1 

The equal-time correlation function ( l52l) is, from ( 156|) and ( 1571) . 

E(r,t)®E(r',t)) = ^ j^^ du coth (^-^^ G,{r,r',co), (58) 

where a factor of 2J\f{u) + 1 has been rewritten as a hyperbolic cotangent (recall (H7|) ). 
Correlation functions for the magnetic field operator are easily found from the electric- 
field expressions through use of (!23|) . with the results 



B^(r, u) ® B(r , u') ) = A7chfioS{uj - u')Af{uj)V x Gi(r, r', u) x V' (59) 

^^"^^ ^B{T,u)®B\r',uj')), (60) 



Af{u) + 1 

B(r,t)®B(r',t)\ = ^ /"~ dw coth f ^^^^ V x Gi(r, r', w) x V' . (61) 
/ TT Jq \2kBT J 

Correlation functions for the reservoir field operators X^^ and Y^^ are found using 
(p^ and fl25|) . Our goal is to calculate the electromagnetic part of the energy density 
and stress tensor, by isolating the electromagnetic part of the expressions ( l33l) and 
fH5]) when they are evaluated in the case of thermal equilibrium. This requires us to 
eliminate the reservoir fields by expressing them in terms of the electromagnetic fields 
and free-field terms that are independent of the coupling to the electromagnetic fields. 

§ Note the distinction between the somewhat confusing terms magnetic media and magnetodielectric 
media. The former refers to media with permanent magnetizations in the absence of appUed fields, the 
latter to media that have a non-trivial magnetic (and electric) response to electromagnetic fields. 
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Equations fl24l) and (1251) perform this separation of the reservoir fields into terms that 
would be present in the absence of any coupling to the electromagnetic fields (the first 
terms on the right-hand sides) and electromagnetic terms (the second terms on the right- 
hand sides). Only terms in correlation functions that depend on the electromagnetic 
fields are required to obtain the electromagnetic part of the energy density and stress 
tensor. We therefore drop terms in correlation functions that are independent of the 
electromagnetic fields. (Further remarks on the non-electromagnetic parts of the energy 
density and stress tensor will be made in the next section.) 

The correlation function in the frequency domain of the reservoir field with 
itself is obtained using IHM . We use a subscript e on the correlation function to denote 
the fact that we include only the terms that depend on the electric field; this yields the 
following correlation function: 

Xt(r,a;')®Xjr',a;'' 




h TTa{r',uj) 
-5(a;-c.) 



UJ — UJ" 



iO^ 



+ 



+ 




2uj 



5{uj 



UJ 



7r« r ,uj\ 



UJ 



1 



1 



Cl{r,uj')(dt{v\uj"] 



+ 



a[T,uj]a[r ,uj] 



X 



Et 



UJ — uj' + iO+ UJ + uj' 
1 



E^l 



+ 



UJ 



uj' + iO+ UJ + uj' 



1 



r, UJ 



+ 



Ce(r',a;"; 



UJ 



UJ" 



iO+ u + w" 



r, UJ 



Efr', 



UJ 



(62) 



The third term on the right-hand side of fl62|) will produce a product of principal values if 
the quantities containing the infinitesimal number O"*" are expanded in terms of principal 
values and delta functions. This product of principal values would need to be treated 
with care depending on the integration variable to which they are referred (as well as 
integrations over uj' and uj" , there will be an integration over uj in the energy density 
and stress). It is safer to rewrite the term in question in a form that does not produce 
a product of principal values; this can be done using the identity 



1 



{uj-uj' + iO+)(w - uj" 



iO- 



UJ' 



UJ" 



2i0- 



uj' - uj" - 2i0 
P 



1 



UJ 



UJ 



UJ 



1 



{uj - uj'){uj - uj") 

5{uj - uj') + 5{uj 



uj' + iO^ 
m5{uj - 



UJ 



UJ" 



iO^ 



UJ 



m5{uj 



UJ 



UJ 



— ivr- 



{uj - uj'){uj - uj") uj'-uj"-2i0+ ' ^^^^ 

From (163|) . the product of the quantities in square brackets in the third term on the 
right-hand side of ( l62l) simplifies to 

1 



UJ — uj' + 



iO+ 



UJ + uj' 



UJ — 
Auj^ 



UJ" 



iO^ 



UJ + uj" _ 
5iuj - uj') + 5iuj 



UJ" 
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where all terms apart from the right-hand side of fl63l) were expanded as principal values 
and delta functions (we do not perform a similar expansion in the last term in fl63|) and 
( 164|) at this stage simply to save space). At the risk of long-windedness, we rewrite ( l62l) 
having inserted fl64|) : 



—6{u-uj) 

2uj uj 




+ 



a[r,u)a[r ,u) 



co-u" - iO+ u + u" 
1 1 



+ 



u — u' + iO+ u + u' 



r, UJ 



w'2)(u;2 



II2\ 



. 6{uj - uj') + 6{uj - uj") 

ivr 

■ 2i0+ 



E(r',a; 



(65) 



The correlation function in the third term on the right-hand side of ( 16511 is given by 
f l56|l : the correlation functions in the first two terms are shown by fH5|) - fj48l) . f l2T|) and 
191) to be 

.1/2 



E\{y,u')C,,{y\uj")) = 27r/iow'W(w' 



he, 



TT 



^i(r,a;") 



G,,{r',r,u")6{u'-u"), (66) 



ei{r',uj') 



1/2 



G*^{r,r\u:')6{io'-io"). (67) 



Now we must focus on the expectation values containing only X^^ that are required 
for the expectation values of the energy density (l33ll and stress tensor (H3|l . In writing 
the expectation values we employ a shortened notation for two limits involving the Green 
bi-tensor: 

A^. .(r,a;) := lim Gij{r,r',u), 



A\.(r,a;) := lim[V x G(r,r',w)x V']ij. 



(68) 
(69) 



The limit r' — )■ r appears because the expectation values of the energy density (l33ll 
and stress tensor (l43l) contain correlation functions evaluated at r' = r. But the Green 
bi-tensor itself diverges when r' = r, so the limit in fl68p and fl69p must be understood 
to be taken only in the final expressions for physical quantities. The zero-point part of 
the Casimir effect requires a regularization to remove the divergent zero-point energy 
that is always present in a homogeneous medium (including vacuum) and that does 
not contribute to the Casimir force. This regularization is implemented at the level 
of the Green bi-tensor in a manner familiar from Lifshitz theory [21 [3l H]. We also 
understand this regularization to be included in fl68l) and fl69l) when they appear below 
in the Casimir energy density and stress tensor. 

For the energy density (!33ll we need dco'(((9X^)2/2 + w^X^ /2); from the general 
relation ( |52l) it therefore follows that we must insert a factor of ibj'bj" + uj'^)/2 into (|65ll . 
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Inserting this factor and substituting (IMj) . (jSZD, dSSD and ^ we obtain, after minor 
simplifications, 



hjoj" + uj')xUv,uj')X^,{v,uj") 



hn (u'u" + 



V^i(r,w)£i(r, 



2lo 



W2 _ uj"^ 



x5iu - uj')6iuj' - uj")A''^,{r,uj') 



, hTT {uj'u" + g;^) ^- 

+ — V w)£i(r, w') 



2uj 



Lip- — Uj'"^ 



— m5{uj — uj') 



h {io'oo" + 



— ITT 



Sjco - tu') + 5(10 - to") 
u' - u" - 2i0+ 



xA/'(u;')'^(w' - w")ImA\.(r, w')- 



(70) 



The right-hand side of (ITOjl consists of a sum of three terms. The first two terms each 
contain a product of three delta functions; these terms containing three delta functions 
can be combined and written as 

'^^^£,{r,u)^^icu) [iA^,,(r,a;') - iA^*^.(r,a;')] S{cu - cu')5icu - cu")5{co' - cu"). (71) 

The third term on the right-hand side of fITU]) has a part containing the sum of two delta 
functions in a numerator; this part gives the contribution 

1 



UJ 



, hix iuj'uj" + 



-ei[r,uj) 



{u' - uj") 
5{uj-uj') 8{uj-uj") 



- + -mdiJ - u") 



+ P- 



(w - cu") {u' - u) 



+ 2i7T6{uj - uj')6{uj - u") 



ei(r, a;)A/'(a;)ImA^,,(r, uj')5{uj - uj')5{uj - uj")5{uj' - u") 



(72) 



which cancels with (17T]) . The first two terms in the sum on the right-hand of (!70!) have 
now been reduced to the parts containing a principal value; consider the first of these, 
namely 



w2 - uj"^ 



M{uj')5{uj - uj')5{uj' - a;")A^ii(r,c^')- (73) 



Because of the delta functions, ( !73|l is restricted to contribute only when the denominator 
in the principal value vanishes; by the definition of a principal value there is no such 
contribution, unless the factor u—u" in the denominator uP—bj"'^ is canceled by an equal 
factor bj — bj" in the numerator. But there is such a factor in the numerator, obtained 
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by Taylor expanding the remaining function of uj (apart from the the factor uj — u" in 
the denominator) around the point u = u"; the contribution of fl73l) is therefore 



h7id[u"ei{r,u")] 



(74) 



The principal value in the second term in the sum on the right-hand of (170|) is dealt with 
in the manner employed for (1731) . Implementing all these simplifications of (170|) . and 
using the symmetry A^^^(r, w) = A^j-{r,uj) that follows from and (ES]), we obtain 



E 



2h'Kd[u'ei{Y,uj')] 



■M{J)b{uj - J)b{J - u;")ReA^,,(r, J) 



(75) 



(76) 



The equality (!76|) is easily verified by tracing back minor changes in the derivation of 
(1751) . We can now write the time-domain expectation value da;((9Xtj)^/2 + a;^X5/2) 
using fl76|) . f l75|) and the general relation fl52|) : 

1 



dw 



1-nc^ 

h 







2kBT 



ReA^/(r,a;) 



dui 







du;' P A my^i(^' ^) coth 



(W2 _ 



no;' 



Use of the identity 



(u;2 - u;'2)2 duj' \u^- u'^ 



ImA^/(r,u;0 (77) 



(7J 



in the last term of f l77|) allows the integration over w to be performed by means of the 
Kramers-Kronig relation ([7]), yielding finally 



do; 



h 



27rc2 



Im 



du' coth 



2kBT I \duj 



;KKr,u;')-l]} A^/(r,a;'). (79) 



Turning now to the the expectation value containing only X^ that appears in 
the expectation value of the stress tensor ( H3|) . we see that this is duj{{dX^Y /2 — 
w^X2/2). Recalling the general relation (152|) . we must therefore insert a factor 
{(jj'uj" — w^)/2 in ( |65ll . This factor, in the first term on the right-hand side of (l65l) . 
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can be written uj'iuj" — w)/2 because of the presence of 5{uj — u'), and so the pole at 
u = u" inside the square brackets is cancelled. Similarly, the pole ai u = u' in the 
second term on the right-hand side of ( l65l) is also cancelled by the factor {u'u" — uj'^)/2. 
In the final term on the right-hand side of (165|) . multiplication by the factor {u'u" —uj'^)/2 
also produces a simplification that follows from 

^ 6{UJ - LO') + 6{UJ - LO") 



-uj'5{uj-uj')+uj"5{uj-u"), 



(80) 



u' - u" - 2i0 

which vanishes when combined with a factor 5{uj' — u") from fl56|) . We substitute (|66|1 
dnZD, ([56]) and ([6]) and find 

1 
2 



[UJ LO 



2h'K 



c^ 



£i(r, uj)M{uj')5{uj - uj')5{uj' - a;")ReA\.(r, J) 



/2 



Y(^(^'^"-^')^-(r,^')^i-(r,u;")) 



(81) 
(82) 



Af{co') + 

The time-domain expectation value dw(((9X(^)^/2 — a;^X^/2) follows from ( IHTl) . ( I82|) 
and the general relation (152|) : the Kramers-Kronig relation ([7]) can be immediately 
applied, with the result 



do; 



2 



h 



27rc2 



Im 



du' coth 



2kBT 



Kr,u;')-l]A^/(r,a;0 



(83) 



Expectation values for the Y^j field analogous to fl79l) and fl83|) are also required. 
In keeping with the discussion earlier in this section, only terms that depend on the 
magnetic-field part of in fl2S]) are included in the calculation and we denote this 
fact by a subscript b on expectation values. The derivations are very similar to those 
described in detail above for the case of the X^^ field; we therefore simply state the 
results: 



do; 



— Im / du' coth 
2% 



hco' 



2kBT J \du 



— {-u'[^{r,u')-l]}]A^:{r,u' 



du 



- — Im / du' coth 
27r 



hco' 



where the definition 



2kBT 
has been employed. 



Kr,a;')-l]A^/(r,a;' 



^5) 
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To calculate the expectation values of the energy density ( 15^ and stress tenor 
we also require the expectation values of duj f3{Y^uj)Y^^iBj and duj a{r,ijj)X^^iEj. 
To find the first of these expectation values, we calculate the frequency-domain 
correlation function of 'Y^^ and B. Again we include only terms that depend on the 
magnetic-field part of 'Y^ in (1251) and denote this fact by a subscript b on the correlation 
function. From f l25|l and f l23|) we find 



Yt(r,a;')®B(r',a;"; 



B 



TX 



'2h. 



-i—\l—S(u-u')(Cl 

V w ' 



V X Wr'tu" 



/5(r,a;) 



2lo 



— 'mb(uj — u') 



V X tHv,uj')®V X E(r',w" 



(86) 



2ujuj'w" 

The second correlation function on the right-hand side of ( l86l) is found from ( l56|) : the 
first correlation function is shown by P5!) - (H8l) . (pTj) and (fT9|) to be 

f?L(r,c.')(V'xE),(r',a;")) 

1/2 



2TTijJ,oUl" 



Ar(u;0(V' X G(r',r,wOx V)jiS{oj' 



(87) 



Insertion of (156|) and ( 1871) in (186!) . with use of ([6]), gives 
tUr,a;')4-(r',a;") 



^ 2 

2vr /ioA/ — 
V w 

27ThflQu' 



TT 

2^0 



1/2 



Ar(w') (V' X G(r', r, u') x V)jiS{io - uj')5{uj' - oo") 



TT 



-ujKi[r, u 



1/2 



5(a;' - u")Ar{u') 



2u 



12 



m5{uj — uj') 



(V X Gi(r,r',a;')x V%-. 



(88) 



It is straightforward to show that the correlation function ^y^j(r, Ci;')-Bj(r', 0;")^ differs 
from fl88l) by a complex conjugation and the replacement of M{uj') by M{uj') + 1; 
we can then compute the equal-time correlation function of and B using these 
frequency- domain results and the general relation (l52l) . The expectation value of interest 
is dcu /3(r, a;)K;i-Bj)B; with use of (l6l) and the Kramer-Kronig relation ([7]), this 
expectation value is found to be 



da;/3(r,u;)F^,(r,t)5,(r,t) 

= Im / duj' [K{r,uj') — 1] coth 



fko' 
2kBT 



A^.(r,a;') 



(89) 



The final expectation value required for our purposes is that of dco a{r,u)Xi^iEj. 
Only terms that depend on the electric-field part of in are included; this fact is 
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denoted by a subscript e on the correlation function. The calculation exactly parallels 
that leading to fl89|) and the result is 



dua{r, u)X^i{r, t)Ej{r, t) 



h 



:Im ^ du' [e{r, u') - 1] coth ^- 



A\-(r,a;') 



TTC- .In - • - \2kBT<- 

This completes the set of expectation values needed to compute the electromagnetic 
part of the energy density and stress tensor in thermal equilibrium. 



6. Casimir energy density 

In the previous section we ignored the free- field parts of the operators and Y^^. 
The rationale for this omission in deriving the Casimir effect is that the stress-energy 
associated with the free-field part of the reservoir represents the absorbed energy due to 
the dissipation of the medium, together with the zero-point energy of the reservoir. The 
dissipated energy is included in the canonical theory of macroscopic electromagnetism, 
with the result that the system is closed and a proper quantization can be performed [10] . 
The description as a closed system is also essential to the existence of a stress-energy- 
momentum tensor, derived in section |3l But Casimir forces are caused by the stress- 
energy of the electromagnetic fields, not by the stress-energy absorbed and dissipated in 
the medium. This is why the correlation functions of the previous section, which include 
only the electromagnetic contribution, determine the Casimir energy density and stress. 

The classical expression f l33|) for the energy density also holds in the quantum 
theory, except that the final term must be written in a hermitian form to give a hermitian 
energy-density operator. We thus have the operator 



+ 



c^ 



doj 



-id A. 



1 



/3 Y^-B + B-Y, 



.(91) 



2 \- 1 — ^ / ' 2^ '' '^/ 2 ^ — "-^ oj/ 2 

The /3-dependent term in fl9T]) has an expectation value that follows immediately from 
( 189|) : the hermitian combination in (19T]) just picks out the real part of (l89l) . which is the 
entire right-hand side. The expectation value of the quadratic terms in the electric and 
magnetic fields in (!9T]) are obtained from (158|) and (!6T!) . When all terms are combined, 
the final expression for the Casimir energy density (p) is 



(p) = — Im / dui coth 
27r Jo 



2kBT 



1 d[uje{r,uj)]^j^^ 



du 



+ 



K r, w — u- 



dK(r, oj) 



du 



A^/(r,a;) 



(92) 



Note that the ^-dependent factor in fl92|) takes the form [d{ufi)/du]/ fi'^ when written 
in terms of p = 1/k. Casimir forces at zero temperature can be found by using fl92|) to 
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calculate the total Casimir energy of a configuration of objects and taking derivatives 
with respect to the parameters specifying their separations and relative orientations. 

The factors in (192|) that depend on the dielectric functions have a familiar form: the 
Brillouin expression [121 E] fo^' the monochromatic electromagnetic energy density of a 
lossless medium contains the same quantities, where e and k are real in that case. The 
result fl92|) is in many ways remarkably simple, considering that it holds for arbitrary 
dispersion (consistent with the Kramers-Kronig relations). Dispersion has a highly 
complicating effect on the electromagnetic energy of general fields in media, even when 
the difficulty of losses can be ignored [15]. Here the losses are compensated because of 
the imposition of thermal equilibrium, and the restriction to thermal (and zero-point) 
fields also has the special effect that dispersion contributes only through the simple 
first-order frequency derivatives in ( !92|) . 

For computational purposes it is more convenient to re-express the frequency 
integral in (!92|) as a sum over imaginary frequencies. Because of the property [1] 

G(r,r',-a;) = G*(r,r',a;) (93) 



of the Green bi-tensor, its real part is even in u while is imaginary part is odd, and the 
same property holds for ^^-ir.uj) and dP^-{r,u) (recall fl68|) and ( 169|) ). Moreover, the 
dielectric functions also have the property ( l93l) [1]. The hyperbolic cotangent in ( l92l) . 
on the other hand, is odd in uj. All this means that the imaginary part of the integral in 
( l92l) is automatically extracted if we modify the integration over u so that it runs from 
— oo to oo and multiply by by — i/2; thus f l92|) can be replaced by 



(P) = -i- 



dw coth 



2kBT J [d 



+ 



1 d[uje{r,uj)]^E 



dnir, uj) 



u- 



du 



A' 



B i. 



(94) 



We can now close the frequency integral in the upper-half complex frequency plane, 
where the dielectric functions and the Green bi-tensor are analytic [TTl |3] . This contour 
integral is given by the sum of the residue contributions from the poles in the hyperbolic 
cotangent term at positive imaginary frequencies u = 2TTikBTn/h, n = 0, 1, 2, . . .. With 
the notation 



2nkBTn 



n 



0,1,2, 



(95) 



for these imaginary frequencies, the Casimir energy density fl94l) then has the form 



(p) = ksTY^ 



n=0 



C 



+ 



1 d[(X'£:(r, o;)] 
'2 



du 



d[ufi{r, uj)] 



doj 



(96) 



where the prime on the summation sign means that the first term in the sum is taken 
with a factor of 1/2 (because the contour passes through the pole at a; = 0), and 
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where k has been replaced by yu(= 1/k). At zero temperature, where only the zero- 
point contribution remains, the sum in f l96|) becomes an integral over positive imaginary 
frequencies: 



The correct expression for the Casimir energy density in media has been a subject of 
conflicting assertions; as pointed out in [16], however, only the form (!97j) gives a Casimir 
force between parallel plates that agrees with the force obtained from the vacuum stress 
tensor between the plates. The correct form of the energy density emerges automatically 
from macroscopic QED, which is moreover a fully quantum treatment of the problem. 

7. Casimir stress tensor 



The quantum stress tensor operator has the same form as the classical expression 
when the latter is written so as to give a hermitian operator: 

^ _ 1 



CTii — -5ij{eQ& + kqW') — EoEiEj — KoBiBj 



+ 



du < 5. 



-{dt±^)' + ^{dtY^)' - ^cu'iXl + Y^) + ( • E + E ■ X, 



(9J 



The Casimir stress tensor is the expectation value of the electromagnetic part of ( l98l) in 
thermal equilibrium. We proceed as in the case of the energy density in the last section. 

The a- and /3-dependent terms in (l98l) have expectation values that follow directly 
from fl90|) and fl89|l . and the expectation values of the terms quadratic in the electric 
and magnetic fields are obtained from ( l58l) and ( l6TI) . The Casimir stress tensor (aij) is 
thereby found to be 



(cTij) = — Im / do; coth 

TT Jo 



2kBT 



c^ 



(99) 



As in the case of the Casimir energy density, a computationally more convenient formula 
for the stress tensor involves a sum over imaginary frequencies. The derivation is as 
described in leading up to fl96|) and the expression is 

'1 

;£^(r,i^„) 



n=0 



25.,A^^(r,ie.)-A\(r,ie. 



+K{r, i^„) 



-<5,,AY(r,ien)-A«^.(r,ia) 



(100) 
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At zero-temperature we obtain from fllOOl) the zero-point Casimir stress as an integral 
over positive imaginary frequencies: 



+/c(r,iO 



-5.,AV(r,iO-A\(r,iO 
^5,,A^^(r,iO-A\.(r,iO 



(101) 



The formula fllOOp (with k = 1) for the Casimir stress tensor in media was obtained 
by Herculean efforts in |3|, and is the most general result of Lifshitz theory. Part of the 
reason why such an enormously complicated formalism was required in [3] is the lack of 
a Hamiltonian and Lagrangian basis for the theory, which also undermines any claims 
that the result applies to quantum electromagnetic fields (see Introduction). As in the 
case of the energy density in the last section, the Casimir stress tensor in media emerges 
in a self-contained manner from macroscopic QED in thermal equilibrium, without the 
need for additional input. 

It is interesting also to compute the expectation value in thermal equilibrium of the 
quantum version of (H4|) . The first term on the left-hand side of (jUj) has of course zero 
expectation value in the stationary situation of thermal equilibrium. The right-hand side 
must be written in hermitian form in the quantum theory, and using aVjO; = Vj(a^)/2 
and /9Vj/3 = Vj(/3^)/2 its expectation value is found in a manner similar to that used to 
obtain (pOj) and (l89l) . This leads to the following result for the divergence of the Casimir 
stress tensor: 



2^^ 



Im 



du coth 



2kBT 



- A^/(r,u;)Vi/t(r, w) 



n=0 



-A 

„2 3 



r, i^„)Vi5(r, i^„) - A^/(r, i^„)ViK(r, i^. 



(102) 
(103) 



As a final remark on the Casimir stress-energy-momentum tensor, we note that 
the energy flux and momentum density in this tensor must be zero, because the 
electromagnetic fields are in thermal equilibrium. It is straightforward to verify, by 
calculations similar to those used to obtain (p2|) and (p9l) . that the electromagnetic 
parts (defined as in section [5l) of the energy flux (IMl) and momentum density (H2l) do 
indeed vanish in thermal equilibrium. 



8. Conclusions 

The Casimir effect has been derived from macroscopic QED [10] by a simple restriction to 
thermal equilibrium. Expressions for the Casimir energy density and stress tensor were 
obtained for arbitrary inhomogeneous magnetodielectrics. As the results are derived 
from a rigorous quantization of electromagnetic fields in dispersive, dissipative media, 
they not subject to the criticisms that have been directed at the standard Lifshitz theory 
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of the Casimir effect. Moreover, the canonical basis of macroscopic QED [TU] means 
that the correct forms of the Casimir energy density and stress tensor in media emerge 
directly from the theory. In Lifshitz theory, by contrast, there is no Hamiltonian or 
Lagrangian, so that detailed mechanical and thermodynamical arguments are required 
to obtain the form of the electromagnetic stress tensor in media. 
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